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ABSTRACT 

We investigate the clustering properties of a complete sample of 10^ star-forming 
galaxies drawn from the data release 4 (DR4) of the Sloan Digital Sky Survey. On 
scales less than 100 kpc , the amplitude of the correlation function exhibits a strong 
dependence on the specific star formation rate of the galaxy. We interpret this as the 
signature of enhanced star formation induced by tidal interactions. We then explore 
how the average star formation rate in a galaxy is enhanced as the projected separation 
Tp between the galaxy and its companions decreases. We find that the enhancement 
depends strongly on r^, but very weakly on the relative luminosity of the companions. 
The enhancement is also stronger in low mass galaxies than in high mass galaxies. In 
order to explore whether a tidal interaction is not only sufficient, but also necessary 
to trigger enhanced star formation in a galaxy, we compute background subtracted 
neighbour counts for the galaxies in our sample. The average number of close neigh- 
bours around galaxies with low to average values of SFR/M* is close to zero. At the 
highest specific star formation rates, however, more than 40% of the galaxies in our 
sample have a companion within a projected radius of 100 kpc. Visual inspection of 
the highest SFR/A/* galaxies without companions reveals that more than 50% of these 
are clear interacting or merging systems. We conclude that tidal interactions are the 
dominant trigger of enhanced star formation in the most strongly star-forming sys- 
tems. Finally, we find clear evidence that tidal interactions not only lead to enhanced 
star formation in galaxies, but also cause structural changes such as an increase in 
concentration. 

Key virords: galaxies: clustering - galaxies: distances and redshifts - large-scale struc- 
ture of Universe - cosmology: theory - dark matter 



1 INTRODUCTION 

It has been known for more than thirty years that 
galaxy interaction s lead to enhanced star formation. 
iToomre fc Toornr3 l| 19721 ) pioneered the use of numerical 
simulations to study the interactions of galaxies and sug- 
gested that gas may be funnelled to the central regions of 
the systems as a result of the strong tidal forces that operate 
during the encounter. This gas is then able to fuel a burst 
of star formation. Since then, there have been many studies, 
both observational and theoretical, that have examined the 
relationship between star formation and galaxy interactions. 

Most e arly observational studi es adopted broad band 
colours (e.g. lLarson fc Tinslevl[r978l ) . Hq equivalent widths. 
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Sanders et al. 19861: Solomon fc S age 19881: iTinnev et ahl 



1990l : lYoung et all Il996l ) as indicators of star formation. 
These studies all demonstrated that galaxy interactions are 
statistically l inked to enhanced rates of star formation (see 
the review of lKeellll99]] : IStruckll 19911 ). 

Recent studies of star formation in interacting galax- 
ies have been based on redshift surveys suc h as the Cen- 
ter f o r Astrophysics red shift survey (CfA2; iBarton et al.l 
I2OO0I : IWoods et allbOOel). the Two D egree Field Redshift 
Survey (2dFGRS; Lambas ct al. 2003'), an d the Sloan Digi- 
tal Sky Survey (SPSS: Nikolic et al. 20o3 : IWoods fc Gelled 
I2OO7I : lEllison et al-llioosl ). These studies have also provided 
observational evidence that star formation is enhanced as 
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a consequence of tidal interactions. Most of these stud- 
ies have also demonstrated that the degree of enhance- 
ment is a strong function of the projected separation be- 
tween the two galaxies as well as their difference in redshift. 
In addition, some studies investig ated how ga laxy prop- 
erties suc h as concentration (iNikolic et al.ll2004l). l uminos- 
ity ratio ll Woods et al.ll2006l: IWoods fc Gelleij|2007l ). steUar 



mass ratio ([E llison et aL| l2008l ). colour, and AGN activity 
l|Woods fc Ge llcr 2007) depend on separation. 

Although most studies have supported the picture that 
interactions induce star formation, t here have been num - 
ber of dissenting papers. For example, iBergvall et al.l (|2003l ) 
analyzed optical/near-IR observations of a sample of 59 
interacting/merging systems and concluded that they do 
not differ very much from isol ated galaxies in te rms of 
their global star formation rates. iBrosch et al.l (|2004l ) found 
that interaction-induced star formation is not significant for 
dwarf galaxies. A more recent study bv lSmith et al.l (|2007l ) 
analyzed Spitzer mid-infrared (MIR) imaging of a sample 
of 35 inter acting galaxy pairs selected from the Arp Atlas 
l| Ard [l96^ ) . They compared the global MIR properties of 
these systems with those of normal spiral galaxies. The MIR 
colors of interacting galaxies were found to be redder than 
normal spirals, implying enhancements to the specific SFRs 
of a factor of ~2. However, in contrast to results from pre- 
vious investigations, they did not find any evidence that the 
enhancement depended on separation. This may be due to 
the small size of their sample an d fact that the ga laxies were 
selected to be tidally disturbed (|Smith et al.ll2007l l. 

On the theoretical side, A^-bo dy simulations that 
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that interactions between galaxies can bring gas from the 
disc to the central regions of the galaxy, lea ding to enhanced 
star f ormation in the bulge. Recently, IDi M atteo et al. 
(|2007l ) investigated star formation in a suite of several hun- 
dred numerical simulations of interacting galaxies with dif- 
ferent gas fractions, bulge-to-disk ratios and orbital param- 
eters. Their work confirmed that galaxy interactions and 
mergers can trigger strong nuclear starbursts. However, the 
authors pointed out that this is not always the case, because 
strong tidal interactions at the first pericenter passage can 
remove a large amount of gas from the galaxy disks, and 
this gas is only partially re-acquired by the galaxies in the 
last phase of the merging event. 

In summary, although it is now well established that 
interactions/mergers between galaxies can enhance star for- 
mation, a number of important questions remain to be an- 
swered: 

• Are interactions not only sufficient but also necessary 
to enhance star formation? 

• Do interactions always trigger enhanced star formation? 

• How does the enhancement in star formation depend on 
parameters such as the separation between the two galaxies 
and their mass ratio? Does the enhancement also depend on 
properties such as stellar mass or galaxy morphology? 

To answer these questions, we adopt three differ- 
ent methods to analyse a sample of ^ 10^ star-forming 
galaxies selected from the Data Release 4 (DR4) of the 



SDSS. First, we compute the cross-correlation between star- 
forming galaxies and a reference sample of galaxies drawn 
from the DR4. In the standard model of structure formation, 
the amplitude of the correlation function on scales larger 
than a few Mpc provides a direct measure of the mass of 
the dark matter haloes that host the galaxies. As we will 
show, the amplitude of the correlation function on scales less 
than ~ 100 kpc can serve as a probe of physical processes 
such as mergers and interactions. We then compute the av- 
erage enhancement in star formation as a function of the 
projected separation between two galaxies and we explore 
how the enhancement depends on galaxy properties such 
as stellar mass and concentration index. Finally, we com- 
pute counts around our galaxies as a function of separation 
and explore how this changes as a function of the specific 
star formation rate SFR/M*. This allows us to investigate 
whether the majority of galaxies with specific star formation 
rates above some critical threshold are experiencing merger- 
induced starbursts. In a separate paper, we explore whether 
AGN activity is also triggered by tidal interactions using the 
same set of analysis techniques. 

Throughout this paper. We assume a cosmological 
model with the density parameter f2o = 0.3 and a cosmo- 
logical constant Ao = 0.7. To avoid the — 51ogiQ/i factor, 
a Hubble constant ft = 1, in units of 100 kms~^Mpc~^, is 
assumed throughout this paper when computing absolute 
magnitudes. 



2 SAMPLES 

2.1 The SDSS Spectroscopic Sample 

The data analyzed in this study are drawn from the Sloan 
Digital Sky Survey (SDSS). The survey goals are to obtain 
photometry of a quarter of the sky and spectra of nearly one 
million objects. Ima ging is obtained in the u, p, r, i, z bands 
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with a special purpose drift scan camera (iGunn et al.lll998f) 
mounted on the SDSS 2.5 meter telescope I Gunn et al.l 
2006) at Apache Point Observatory. The imaging dat a are 
photometrically (iHogg et al.ll200ll : Tfucker et al.1 [20061 ) and 
astrometrically (|Pier et al.l bOOsf ) calibrated, and used to 
select stars, galaxies, and quasars for follow-up fibre spec- 
troscopy. Spectroscopic fibres are assigned to objects on the 
sky using an efficient tiling algori thm designed to optimize 
completeness l|Blanton et al.ll2003l). The de tails of the survey 
strategy can be found in lYork et al.l (|2000l ) and an overview 
of the data pipelines and products is provid ed in the Early 
Data Release paper (jStoughton et al.ll2002l ). More details 
on th e photometric pipeline can be found in iLupton et al.l 
l|200ll ). 

Our parent sample for this study is composed of 397,344 
objects which have been spectroscopically confirmed as 
galaxies a nd have da ta publicly available in the SDSS Data 
Release 4 (jAdelman- McCarthy et al. 2006). These galaxies 
are part of the SDSS ' main' galaxy sampl e used for large 
scale structure studies l|Strauss et al.ll2002l ) and have Pet- 
rosian r magnitudes in the range 14.5 < r < 17.77 after 
correction for fore ground gala c tic ex tinction using the red- 
dening maps of iSchlegel et al.l l|l998t ). Their redshift distri- 
bution extends from ~ 0.005 to 0.30, with a median z of 
0.10. 
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The SDSS spectra are obtained with two 320-fibre spec- 
trographs mounted on the SDSS 2.5-meter telescope. Fibers 
3 arcsec in diameter are manually plugged into custom- 
drilled aluminum plates mounted at the focal plane of the 
telescope. The spectra are exposed for 45 minutes or un- 
til a fiducial signal-to-noise (S/N) is reached. The median 
S/N per pixel for galaxies in the main sample is ~ 14. The 
spectra are processed by an automated pipeline, which flux 
and wavelength calibrates the data from 3800 to 9200 A. 
The instrumental resolution is R = \/S\ = 1850 - 2200 
(FWHM~ 2.4 A at 5000 A). 



2.2 Star- forming galaxies 

Our sample of star-forming galaxies is drawn from the 
DR4 spectroscopic samp le using the criteria described in 
iBrinchmann et al] (|2004[ l. In order for a galaxy to be se- 
curely classified as star-forming, we require that the four 
emission lines [OIII], H/3, Ha and [Nil] all be detected with 
signal-to-noise greater than 3 and that the ratios [OIII] /H/3 
and [NII]/ Hq have values that place them within the re- 
gion of the lBaldwin et all (|l98ll , BPT) diagram occupied by 
galaxies in which the primary source of ionizing photons is 
from HII regions rather than an AGN. We refer to this sam- 
ple as the high S/N star-forming class. In certain cases, we 
supplement the sample with the low S/N star- forming class 
defined by Brinchmann et al. These are the galaxies that 
are left over after all the AGN and high S/N star- forming 
galaxies have been removed, and they have S/N > 2 in Ha. 
Star formation rates can still be estimated from their emis- 
sion line strengths, but the errors on these estimates will be 

significantly larger than for t he high S/N sample . 

The reader is referred to IBrinchmann et al] (|2004l ) for 
a detailed description of how star formation rates are de- 
rived for the various samples. We will be making use of the 
specific star formation rate SFR/M^ estimated within the 
3 arsecond SDSS fibre aperture. These star formation rates 
are more accurate than the total star formation rates derived 
by Brinchmann et al, because they depend only on the emis- 
sion line fluxes measured from the spectra and they do not 
involve any uncertain colour corrections. The disadvantage 
of the fibre-based specific star formation rates is that they 
are only sensitive to the emission from the inner region of 
the galaxy, which includes one third of the total light on 
average. 



2.3 Reference Samples 

We work with two different reference samples: (i) a spec- 
troscopic reference sample, which is used to compute the 
projected cross-correlation function Wpifp) between star- 
forming galaxies and reference galaxies, and (ii) a photo- 
metric reference sample, which is used to calculate counts 
of close neighbours around star-forming galaxies. We use 
the New York University Value Added Galaxy Catalogue 
(NYU-VAGC) to construct the reference samples. The orig- 
inal NYU-VAGC is a catalo gue of local galaxies ( mostly be- 
low z ~ 0.3) constructed by Blanton et al.l (|2005l ) based on 
the SDSS DR2. Here, we use a new version of the NYU- 
VAGC (Sample dr4), which is based on SDSS DR4. The 
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NYU-VAGC is described in detail in lBlanton et al.1 (120051 ). 



Figure 1. Projected redshift-space 2-point cross-correlation func- 
tion tOp(rp) between star-forming galaxies and the reference 
galaxy sample. Different lines correspond to star-forming galax- 
ies with different specific star formation rates. See the text for a 
detailed description. 



The reference samples are exactly the same as used in 
iLi et al. I (l2006b). In short, the spectroscopic reference sam- 
ple is constructed by selecting from Sample dr4 all galaxies 
with 14.5 < r < 17.6 that are identified as galaxies from 
the Main sample, in the redshift range 0.01 2 ^ 0.3, and 
with absolute magnitudes —23 < Mo i^ < —17. The spectro- 
scopic reference sample contains 292,782 galaxies. The pho- 
tometric reference sample is also constructed from Sample 
dr4 by selecting all galaxies with 14.5 < r < 19. The result- 
ing sample includes 1,065,183 galaxies. In certain cases, we 
will work with photometric reference samples with a range 
of differing limiting magnitudes. 



3 CROSS-CORRELATION FUNCTIONS 

Our methodology for computing correlation funct ions has 
been described in detail in our previous papers (|Li et al.l 
l2006al lb[). We present here a brief description and the reader 
is referred to the earlier papers for details. Random sam- 
ples are constructed that have the same selection function 
as the reference sample. The redshift-space two-point cross- 
correlation function (2PCCF) ^(rp,7r) between star-forming 
galaxies and the reference s ample is t hen ca lculated using 
the estimator presented in iLi et al.l l|2006b| j. Finally, the 
redshift-space projected 2PCCF Wp{rp) is estimated by in- 
tegrating (,{rp,TT) along the line-of-sight direction n with |7r| 
ranging from to 40 Mpc. We have also corrected care- 
fully for the effect of fibre collisi ons and a description and 
tests of the method are given in iLi et al. 1 (2006b,i . The er- 
rors on the clustering measureme nts are estimated us ing the 
bootstrap resampling technique (|Barrow et al.|[l984l ). 
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Figure 2. Similar to Figure 1, but in different intervals of stellar mass as indicated at the top of the figure. The symbols are the same 
as in Figure [T] except that a power law corresponding to ^(r) = {r /5h~^ Mpc)~^-^ is additionally plotted in each panel as a long-dashed 
line. 



9.5<logi(,(M/Mg)<10.0 10.0<log,o(M/Mo)<10.5 10.5<log,(,(M/Mg)<11.0 11.0<log,o(M/Mg)<11.5 




0.01 0.1 1 10 0.01 0.1 1 10 0.01 0.1 1 10 0.01 0.1 1 10 
Tp [h-'Mpc] Tp [h-'Mpc] Tp [h-'Mpc] [h-'Mpc] 



Figure 3. Same as Figure[2l except that the three SFR/Mt, samples in each panel are matched in concentration. 



We first compute ■Wp{rp) for our sample of high S/N 
star-forming galaxies from the SDSS DR4. In order to study 
how this depends on star formation rate (SFR), we rank all 
the high S/N star- forming galaxies according to the values 
of their specific star formation rates (SSFR), SFR/M,, and 
define subsamples of 'high SSFR' and 'low SSFR' galaxies as 
those contained within the upper and lower 25th percentiles 
of the distribution of this quantity. The results are shown in 
Figure [T] The dashed (dotted) line corresponds to the high 
(low) SSFR subsample, while the solid line shows the result 
for the sample as a whole. 

Figure [1] shows that galaxies with higher SFR/M^ have 
stronger clustering on scales smaller than 0.1 Mpc and the 
effect becomes stronger at smal ler projected separations. As 
pointed out bv lLi et al] (|2006bl) . the clustering amplitude of 
galaxies depends on a variety of galaxy properties, includ- 
ing stellar mass and galaxy structure. If we wish to isolate 
the effect of the specific star formation rate, it is impor- 
tant that we make sure that the galaxy samples that we 
study are closely matched in terms of other properties, so 
that the effect on the star formation rate can be isolated. 



We have thus divided all the high S/N star-forming galax- 
ies into four subsamples according to logjQ(M,/M0). For 
each subsample we repeat the clustering analysis as above. 
The results are shown in Figure [2] The four panels cor- 
respond to different intervals of logj^Q(M»/M0). To guide 
the eye, a power law corresponding to a real-space 2PCF of 
5(r) = (r/5/i~^Mpc)~^ * is plotted as a long-dashed line in 
each panel. We see that the amplitude of Wpij-p) increases for 
galaxies with larger stellar masses. This is consistent with 
our previous findings about the mass dependence of galaxy 
clustering. We also see that the difference in clustering be- 
tween galaxies with high and low SFR/M, on scales smaller 
than 0.1 Mpc is most pronounced in the lowest stellar mass 
interval. Next, in each of the four M* intervals, we match 
the three SFR/Ma, samples in concentration parameter C 
by requiring that the distribution of C is exactly the same 
as in each of these samples. The Wp{rp) measurements for 
such matched samples are shown in Figure |3] The results 
are very similar to those shown in the previous figure. 

We conclude that the small scale clustering dependences 
shown in Figure [T] are genuinely related to the differing spe- 
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Figure 4. The projected 2PCCF Wp{rp) normalized by the power law corresponding to a real-space 2PCF of ^(r) = {r/5h^^ AIpc)^^'^ , 
as measured at different physical scales and as a function of SFR/Mr. Dashed lines are for high S/N star-forming galaxies only, while 
solid lines show results for the sample including both high and low S/N star-forming galaxies. 



cific star formation rates of the galaxies in the diflterent sam- 
ples. Galaxies with the highest specific star formation rates 
apparently have an excess of companions on scales less than 
100 kpc when compared to the average star-forming galaxy. 
The fact that the increase in clustering occurs only on very 
small scales suggests that the excess star formation is being 
triggered by tidal interactions with these companions. An- 
other intriguing result shown in these figures is that galaxies 
with low SFR/Mt are less clustered on small scales. This 



suggests that there might be a continuous trend linking av- 
erage number of close neighbours and SFR/M,. 

To investigate this in more detail, we calculate how the 
clustering amplitude depends on SFR/M, at a variety of dif- 
ferent physical scales. The results are shown in Figure 3] (red 
dashed lines) . One problem with the high S /N star-forming 
sample is that it does not extend to SFR/Mt values much 
below ~ —10.5. To extend our analysis to lower values, we 
in clude the sample of low S/ N star- forming galaxies defined 
bv lBrinchmann et al.l l|2004h . As discussed in section 2.2, the 
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star formation rates in these galaxies are estimated from the 
Ha hne luminosity, but the dust correction is quite uncertain 
because H/3 is not usually detected with high S/N. Results 
where the low S/N star- forming galaxies have been included 
are plotted as black solid lines in Figure 2] As can be seen, 
logSFR/Mt extends down to values around ~ — 11 for this 
sample. 

On scales larger than 100 kpc, there is very little depen- 
dence of clustering amplitude on specific star formation rate 
for \ogSFR/M, > -10. At lower values of SFR/M^, the 
clustering amplitude increases. This is a manifestation of the 
strong relation between star formation and local density or 
environment. It is well-known that galaxies located in dense, 
massive structures such as clusters have l ower specific star 
form ation rates than "field" galaxies (e.g. iKauffmann et ahl 
|2004 ). It is currently accepted that after a galaxy is accreted 
onto a larger structure, such as a group or cluster, its star 
formation rate will decline, either because its gas is removed 
by processes such as ram-pressure stripping, or simply be- 
cause no further gas accretion takes place and the galaxy 
runs out of the fuel to make new stars. 

On scales less than 100 kpc, the dependence of cluster- 
ing amplitude on specific star formation rate is more com- 
plicated. At values of log SFR/M^ less than -10, we see the 
same increase in clustering amplitude that we saw on larger 
scales. This may appe ar somewhat surprisi ng at first. In a 
recent paper, however, iBarton et al.l (|2007l ) use cosmologi- 
cal simulations to show that a substantial fraction of galax- 
ies selected as "close pairs" from surveys such as SDSS or 
2DF, do in fact reside in very massive dark matter halos. 
Based on this work, we conjecture that the rise in clustering 
amplitude seen at all separations at low values of SFR/M, 
rate is the result of star formation shutting down in galaxy 
groups and clusters. At values of \ogSFR/Mt greater than 
-10 and at separations less than ~ 50 kpc, the clustering 
amplitude shows a strong and continuous increase towards 
larger values of SFR/hh. This is a clear signal that mergers 
or interactions play an important role in triggering enhanced 
star formation in galaxies. 



4 STAR FORMATION ENHANCEMENT 
FUNCTIONS 

In this section, we probe the relationship between star for- 
mation and galaxy interactions by quantifying the enhance- 
ment in star formation as a function of the projected separa- 
tion between two galaxies. We also study how the enhance- 
ment depends on the physical properties of the main galaxy. 
We compute how the average value of SFR/M^ changes as 
a function of the projected distance to the neighbours. From 
now on, we restrict our attention to the sample of high S/N 
star-forming galaxies. These lie in low density environments 
where processes such as ram-pressure stripping, gas starva- 
tion etc should play a much less important role (see Figure 
4). 

The neighbours of a galaxy are identified using the pho- 
tometric reference sample. The advantage of using the pho- 
tometric sample is that the result is not affected by incom- 
pleteness (e.g. the effect of fibre collisions). However, the dis- 
advantage is that some fraction of the close neighbours will 
not be true nearby systems, but rather chance projections of 



foreground and background galaxies that lie along the line- 
of-sight. We correct for this as follows: We count the number 
of companions in the photometric reference sample at a pro- 
jected physical distance for each galaxy with a high S/N 
measure of the specific star formation rate \og{SFR/M,). 
We also generate 10 random samples that have the same 
geometry as the photometric reference sample by random- 
izing the sky position of the photometric objects and keep- 
ing all the other quantities (e.g. the magnitudes) fixed. We 
use these random catalogues to estimate the mean num- 
ber of projected companions expected at random around 
each galaxy. The true number of companions at separation 
rp is given by the difference between the observed and the 
projected random companion count. We then calculate a 
weighted average specific star formation rate at projected 
distance rp by weighting each galaxy by its true compan- 
ion number. The enhancement in \og{SFR/AI,), Ex{rp), is 
defined as the difference between the weighted average and 
the unweighted one. This can be written as 

77 / N Y.'^ ^AnoA{rp) - np,i{rp)] Xi 

^x(rp) = ^jvr — ~ r~r-, ]y — ' ^ ' 



J2i [^o,i{rp) - np,,{rp)] 



where Xi = \og{SFRi/Mt,i) is the specific star formation 
rate of the i'th galaxy, and no,i and the observed and 

projected random companion counts as described above. 

We first consider all high S /N star-forming galaxies with 
r-band apparent magnitude in the range 14.5 < r < 17.6. 
To begin, we restrict the photometric reference sample to 
galaxies with r-band magnitudes r < 19.0. In order to ensure 
that we are finding similar neighbours at all redshifts, we 
only consider neighbouring galaxies that are brighter than 
rsFG + 1.4 mag. The result is shown in Figure[5] The errors 
are estimated using Bootstrap resampling techniques. The 
dashed lines indicate the variance between 10 samples in 
which we randomize the sky positions of the star-forming 
galaxies. 

The top left panel of Figure [S] gives the result for the 
sample as a whole. On scales larger than a Mpc or so Ex is 
constant at a slightly but significantly negative value. This 
is because the average in equation 1 is pair-weighted, and 
galaxies in massive halos have lower specific star formation 
rates than average, but more "companions" at large rp than 
average as a result of large-scale bias effects. On scales below 
about 100 kpc. Ex increases sharply and reaches values cor- 
responding to a factor of about two at projected separations 
less than 20 kpc. In the top right panel of Figure [S] we plot 
results for galaxies divided into two different ranges in stellar 
mass. These results show that there is a strong dependence 
of star formation enhancement on galaxy mass, in that star 
formation in small galaxies is more strongly enhanced at a 
given projected separation. 

In the bottom panels, we scale the projected separation 
rp by the physical size of the galaxy. We use iigo, the ra- 
dius containing 90% of the total r-band light, to calculate a 
scaled projected separation and recompute the enhancement 
function as a function of this scaled quantity. One can see 
that the results are quite similar. Star formation is enhanced 
at separations less than ~ 10 times the optical radius of the 
galaxy and the effect is stronger for lower mass systems. 

We now investigate the importance of the relative mass 
of the companion galaxy in determining the degree to which 
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Figure 5. Star formation enhancement as a function of the projected separation Vp (top panels) and as a function of the scaled separation 
^'p/R'M (bottom panels), for all the high S/N star forming galaxies (left panels) and for galaxies in different stellar mass ranges (right 
panels). All the errors are estimated using the Bootstrap resampling technique. The dashed lines in each panel indicate the variance 
between 10 realizations in which the sky positions of the star-forming galaxies are randomized. See the text for details. 



star formation is enhanced in the primary galaxy by an- 
alyzing galaxy samples with different limiting magnitudes. 
We first keep the magnitude limit of star-forming sample 
constant at rsFG = 17.6, and explore what happens if we 
change the limiting magnitude of the reference sample. The 
results are plotted as circles for rpho < 18.5 and as trian- 
gles for Tpho < 19.5 in the left-hand panels of Figure [6] The 
result shown in the previous figure is plotted as squares. 
Next, we fix the limiting magnitude of the reference sam- 
ple at rpho = 19.5, but decrease the magnitude limit of 



the star-forming sample. The maximum allowed difference 
in magnitude between the star-forming galaxy and its com- 
panion is also increased accordingly. Results are shown for 
rsFG < 16.5, 17.0 and 17.6. in the right-hand panels of Fig- 
ure [S] 

We see from Figure |6] that the star formation enhance- 
ment depends very little on the mass ratio between the star- 
forming galaxy and its companion. There are small changes 
in the expected direction (i.e. there is slightly less enhance- 
ment for companions with lower relative mass), but to first 
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Figure 6. Star formation enliancement as a function of projected separation rp (top panels) and as a function of the scaled separation 
TpjHiM (bottom panels). In the left-hand panels, different symbols connected by solid lines correspond to reference samples with different 
limiting magnitudes (as indicated), while the magnitude of star-forming samples is kept constant at rgpQ = 17.6. In the right-hand 
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All the errors are estimated using the Bootstrap resampling technique. The dashed lines in each panel indicate the variance between 10 
realizations in which the sky positions of the star-forming galaxies are randomized. See the text for details. 



order the enhancement function remains remarkably con- 
stant for different mass ratios. In Figure [T] we divide the 
star-forming sample into two different stellar mass intervals 
and explore if our results change. We find that the enhance- 
ment function has very little dependence on the mass ratio 
of the companion for both low mass and high mass star- 
forming galaxies. 

Finally, we investigate the enhancement function for 



galaxies with different structural properties. We divide all 
the high S/N star- forming galaxies into different intervals 
of concentration parameter C and repeat the analysis de- 
scribed above for each of these subsamples. In Figure |8] we 
plot the results for star-forming galaxies with tsfg < 17.6 
and for reference galaxies with Vp^o < 19.5. We see that the 
star formation enhancement does depend on C, in that the 
galaxies with larger C values are more strongly enhanced. 
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Figure 7. The same as the bottom- left panel of Figure |6] but for the low mass (the left-hand panel) and the high mass (the right-hand 
panel) subsamples separately. To guide the eye, the result for the whole sample in case of rsFG < 17.6 and Tp^o < 19.0 is plotted as 
solid black circles in every panel. 




Figure 8. Similar to the previous plot but for star-forming galaxies with different concentration indices, as indicated above each panel, 
and for the case of rsFG < 17.6 and rp^o < 19.5 only. The blue triangles are for the concentration subsamples and the black solid circles 
are for the whole sample. 



One possible explanation for this effect is that interaction- 
induced starbursts occur when gas flows into the core of a 
galaxy, causing it to bec ome more centrally concentrated 
(I Sanders fc MirabelllToQ^ . 

In Figure[9l we investigate how the concentration index 
of a star-forming galaxy depends on its location in the plane 
of specific star formation rate versus stellar mass. We see 
that at fixed stellar mass, the average concentration index 
is highest for galaxies that are currently experiencing both 
higher-than-average and lower-than-average rates of star for- 
mation. One interpretation of this plot is that tidal interac- 
tions cause gas to fiow from the disk to the nucleus and this 
triggers a starburst at the centre of the galaxy. The forma- 
tion of new stars in the central regions causes the concen- 
tration index to increase. The starburst is then followed by 



a period of relative quiescence, which lasts until the galaxy 
is able to accrete more gas into its disk. Formation of stars 
in the disk brings the galaxy back into the "central plane" 
occupied by galaxies with logSFR/M* ~ —9.5 in Figure |9l 
In the bottom-right corner of the plot, both stellar mass and 
concentration are high, but the specific star formation rate 
is low. This is the regime of early-type galaxies. 



5 CLOSE NEIGHBOUR COUNTS 

In this section, we investigate whether tidal interactions are 
not only a sufficient, but also a necessary condition for a 
galaxy to experience enhanced star formation. We count the 
number of galaxies in the photometric sample in the vicinity 
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Figure 10. Average counts of galaxies in the photometric sample (panels from left to right: rj^^ < 18, 19, and 20) within a given 
projected radius Rp from the star-forming galaxies. Different symbols are for star-forming galaxies in different intervals of specific star 
formation rate, as indicated. 
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Figure 9. Distribution of both high and low S/N star-forming 
galaxies in the plane of stellar mass versus specific star formation 
rate, coloured by concentration index iJ90/R50 measured in the 
5;-band. The color coding of iJ90/i?50 is shown in the bar at the 
right-hand. 



of the star-forming galaxies and make a statistical correction 
for the effect of chance projections by subtracting the aver- 
age count around randomly placed galaxies. 

In Figure [10] we plot the average correlated neighbour 
count (i.e. after statistical correction for uncorrelated pro- 
jected neighbours) within a given value of the projected ra- 
dius Rp. Results are shown for high S/N star-forming galax- 
ies in different intervals of specific star formation rate. We 
have trimmed each subsample so that they each have the 
same distribution in redshift and in stellar mass M« . Panels 
from left to right correspond to photometric reference sam- 
ples that are limited at r = 18.0, 19.0 and 20.0. The star- 
forming sample always has a limiting magnitude of 17.6. 

Figure[TD]shows that the counts around the star-forming 
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Figure 11. Same as the right-hand panel of Figure [TOl but for 
289 galaxies that have the highest specific star formation rates 
{logiQ{SFR/Mf) > —8.8). Results are shown only for scales be- 
low lOOkpc. 



galaxies with different specific star formation rates match 
well on large scales. On scales smaller than ~ 100 kpc, there 
are strong trends in the number of neighbours as a function 
of SFR/M*; galaxies with higher star formation rates are 
more likely to have a near neighbour. 

It is interesting that the average number of close 
neighbours around galaxies with low-to-average values of 
SFR/Mt is close to zero on scales less than 20-30 kpc. On 
scales less than 100 kpc, only around 3% of the galaxies 
in the lowest SFR/A/, bin have a companion. This implies 
that tidal interactions that do not result in enhanced star 
formation are a rare occurrence. This is consistent with the 
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Figure 12. SDSS optical image for 20 star-forming galaxies that are included in the subsample of highest specific star formation 
rate(log]^Q(S'Fi?/Mt) > —8.8). Images are shown for 10 galaxies that are classified as "mergers" (top panels) and for 10 galaxies that 
are classified as "no mergers" (bottom panels). 



findings of lOi Matteo et al.l (|2007h who find that about 85% 
of their simulated sample of interacting galaxies show an en- 
hancement in star formation by a factor > 2. As the specific 
star formation rate increases, the average number of close 
neighbours also rises. Fifteen percent of galaxies in our high- 
est SFR/M» bin have a companion within 20-30 kpc and 
this rises to a value close to 30% if we consider compan- 
ions within 100 kpc from the primary galaxy. In Figure [TTl 
we show the result for 289 galaxies with the very highest 
specific star formation rates in our sample. The fraction of 
galaxies that have a companion within 100 kpc increases to 
values close to 40%. However, it is still true that not every 
star-bursting galaxy in our sample has a close companion. 

How can we explain those galaxies with very high 



SFR/Mf but no close neighbours? We have visually exam- 
ined the SDSS r-band images of 160 star-forming galaxies 
that are included in the highest SFR/M, subsample but 
have no companions within 50 kpc. We classified the sys- 
tems according to whether or not they show clear signs of 
mergers or interactions, including double nuclei and tidal 
tails. Three of us (CL, GK and Roderik Overzier) did this 
independently to make sure that we obtained the same an- 
swer. We find that more than half of such galaxies show clear 
evidence of recent mergers. In Figure [T^ we show some ex- 
amples of our classifications. We thus conclude that, at least 
for galaxies with very high SFR/AI,, interactions or mergers 
are the dominant mechanism for triggering and enhancing 
their star formation. 
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Figure 13. Left: Ope circles show the fraction of star-forming galaxies that have at least one companion in the photometric sample 
within a given projected radius Rp. The solid line shows the result after the number of companions is corrected using random samples. 
Right: the fraction of "isolated" (red) and "non-isolated" (blue) paired galaxies, classified according to whether they have companions in 
the spectroscopic reference sample with the projected separation rp in the range 100 /i~^kpc < rp < Rp and velocity difference smaller 
than 500 km s~^. The parent sample of paired galaxies is selected from all the high S/N star-forming galaxies by requring that a paired 
galaxy has at least one companion within a projected radius of 50 /i~^kpc in the photometric reference sample limited at Vpho = 19.0. 
Solid/dashed lines (open circles/squares) show results obtained with the corrected (observed) companion numbers. 
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Figure 14. Distribution of specific star formation rate for the "isolated" (solid) and "non- isolated" (dashed) paired galaxies. In the 
right-hand panel, the two samples are matched closly in stellar mass. 
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6 EFFECT OF RICH ENVIRONMENTS 

Fig.[5](§|4)) shows that the enhancement in log{SFR/M,) on 
large scales is not zero. Rather, it is constant at a small but 
significantly negative value. We have attributed this large- 
scale bias to the well-known anti-correlation of the star for- 
mation rates in galaxies with the richness of their local envi- 
ronment. The question then arises whether the star forma- 
tion enhancements that we compute may not be a true reflec- 
tion of the effect of galajcy-galaxy interactions, but may be 
biased because some close pairs are not real interacting sys- 
tems, but are associated with group/cluster environments. 

In order to address this problem, we have selected a 
sample of 9052 paired galaxies from all the high S/N star- 
forming galaxies. A galaxy is defined to be paired if it has at 
least one companion within a projected radius of 50 h~^kpc 
in the photometric reference sample limited at Vpho ~ 19.0. 
We use the the photometric sample rather than the sepectro- 
scopic one to select pairs, so that the resulting sample is not 
biased to pairs with approximately equal masses. The max- 
imum pair separation is chosen to yield a sample in which 
the contamination by chance projections is negligible. This is 
demonstrated in the left-hand panel of Fig. [TS] The fraction 
of star-forming galaxies that have at least one companion 
in the photometric sample within a given projected radius 
Rp is plotted as a function of Rp as open circles, and is 
compared to the result (the solid line) after the number of 
companions is corrected with the help of random samples 
(see §|4]for details). As can be seen, at separations smaller 
than 50 h'^kpc, the correction is negligible. 

Next, We classify the paired galaxies as either "isolated" 
or "non-isolated" according to whether they have compan- 
ions in the spectroscopic reference sample with projected 
separation rp in the range 100 /i~^kpc < rp < Rp and veloc- 
ity difference smaller than 500 km s^^. A galaxy is isolated 
if this ring-like region is completely empty. In contrast, a 
non-isolated galaxy is required to have at least one com- 
panion in the annulus that is brighter than that galaxy in 
the r-band. The right-hand panel of Fig. [T3] shows how the 
fraction of these two types of paired galaxies changes as the 
maximum separation Rp increases. We see that the sam- 
ple is always dominated by isolated pairs. However, even 
a small fraction of non-isolated galaxies may still change 
the enhancement function significantly, because the quanti- 
ties defined in equation 1 are pair-weighted and galaxies in 
richer environment have more companions. In the left-hand 
panel of Fig. [14] we compare the distribution of specific star 
formation rates for the isolated and the non-isolated sam- 
ples. The right-hand panel shows the results after the two 
samples are matched in stellar mass with each other (This is 
important because m assive galaxies a re more clustered than 
less massive galaxies (|Li et al.ll2006al )). As can be seen, the 
effect of a rich environment is quite small. The average spe- 
cific star formation rate of the isolated sample differs from 
that of the non-isolated sample by < 0.03 dex only. This 
is consistent with the findings of iBalogh et al.l (|2004| ) who 
showed that the main effect of a dense environment on large 
scales is that the fraction of emission-line galaxies decreases, 
but that the distribution of equivalent widths of Ha among 
the star- forming population remains virtually constant. Here 
we show that the same thing holds for galaxies with close 
pairs. We thus conclude that the large-scale environment 



does not affect the measured star formation enhancements 
in our sample of high S/N star- forming galaxies on scales 
below 100 h-^kpc. 



7 SUMMARY AND DISCUSSION 

We find that the clustering amplitude of high S/N star- 
forming galaxies depends strongly on the specific star forma- 
tion rate on scales less than 100 kpc. The clustering ampli- 
tude increases smoothly as a function of SFR/M* and the 
increase in amplitude is largest at the smallest projected 
separations. We interpret this behaviour as the signature of 
tidal interactions, which lead to infiow of gas and an en- 
hancement in star formation in the two interacting galaxies. 
At low values of SFR/M*, the clustering amplitude again 
increases. The increase occurs on all scales and probably re- 
flects the fact that star formation in galaxies switches ofi' 
after they are accreted onto larger structures such as groups 
and clusters. 

We have explored how the average star formation rates 
of galaxies are enhanced as a function of the projected sep- 
aration of their companions. The enhancement is a strong 
function of separation, increasing from zero at rp > lOOkpc 
to factors of 1.5-4 at rp — 20 kpc. We find that the enhance- 
ment at given separation is stronger for lower mass galaxies. 
Remarkably, we find that the enhancement has almost no 
dependence on the relative luminosity of the companions. 

The tidal force between two objects is expected to scale 
as and m/M, where d is the separation and m/M is 
the mass ratio between the two objects, so it is perhaps not 
surprising that we see stronger star formation enhancement 
as a function of separation than as a function of mass ratio. 
Nevertheless, in this study, we find that the effect of a com- 
panion that is 3 magnitudes fainter than the primary galaxy 
is very similar to that of a companion that is only a factor 
of 3 less luminous than the primary. This is quite startling 
and is worthy of further investigation. 

In order to explore whether tidal interactions are not 
only sufficient, but also a necessary condition for enhanced 
star formation in a galaxy, we have computed background 
subtracted neighbour counts around the galaxies in our sam- 
ple. We find that the average number of galaxies around 
galaxies with low values of SFR/M* is very small. At the 
very highest specific star formation rates, more than 40% 
of the galaxies in our sample have a companion within a 
projected radius of 100 kpc. Visual inspection of the high 
SFR/Af* galaxies without companions reveals that more 
than 50% of these are clear interacting or merging systems. 
We thus conclude that tidal interactions are the primary 
mechanism for inducing the highest rates of star formation 
in galaxies in the local Universe. 

Finally, we find clear evidence that tidal interactions not 
only lead to enhanced star formation, but also cause struc- 
tural changes in galaxies. Many of the most strongly star- 
forming galaxies in our sample have concentration indices 
similar to those of normal early- type galaxies. We note that 
the concentration index is measured for the r-band light and 
not for the stellar mass, so more careful analysis is needed 
before one can definitely conclude that the interactions will 
result in the formation of a galaxy with a high bulge mass 
fraction. Nevertheless, we conclude that our results are in 
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general accord with the t heoretical picture first laid out by 
iToomre &: Toornr3 l| 19721 ). which showed that galaxy inter- 
actions can lead to the growth of bulges and spheroids in 
galaxies. 
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